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We construct a mean-field theory for itinerant ferromagnetism coexisting witii a non-unitary superconducting 
state, where only the majority-spin band is gapped and contains fine nodes, while the minority-spin band is 
gapless at the Fermi level. Our study is motivated by recent experimental results indicating that this may be the 
physical situation realized in the heavy-fermion compound UGe2 . We investigate the stability of the mean-field 
solution of the magnetic and superconducting order parameters. Also, we provide theoretical predictions for ex- 
perimentally measurable properties of such a non-unitary superconductor: the specific heat capacity, the Knight 
shift, and the tunneling conductance spectra. Our study should be useful for direct comparison with experimental 
results and also for further predictions of the physics that may be expected in ferromagnetic superconductors. 
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I. INTRODUCTION 

The interplay between ferromagnetic (FM) and supercon- 
ducting (SC) long range order microscopically coexisting in 
the same material has attracted much interest during the last 
decade due to the discovery of superconductivity in ferro- 
magnetic metals, UGe2, URhGe, UCoGe, and possibly ZrZn2 
(see, however, Ref. j^ .^A^A^ . One possible route of investi- 
gation of such systems was adopted in early works^-^ '^, which 
assumed a conventional s-wave superconducting condensate 
residing in a ferromagnetic background caused by localized 
spins or aligned magnetic impurities. It was shown that be- 
low a critical value of the magnetic coupling, comparable to 
the superconducting gap A itself, superconductivity and mag- 
netism were able to coexist. It was also suggested that a finite 
momentum pairing state, known as FFLO phase can appear 
in the presence of external magnetic field or intrinsic ferro- 
magnetic order, and could thereby permit a larger threshold of 
the spin exchange energy to coexist with superconductivity. 

On the other hand, it has been known since the early days of 
research on -^He that alternative superconducting states, other 
than s-wave, can be favoured in a ferromagnetic background. 
The early theories of an equivalent phenomenon to occur in 
the solid-state have been formulated in the early 1980s", de- 
spite the absence of any experimental example of a ferromag- 
netic superconductor at the time'^. With the discovery of su- 
perconductivity in UGe2 and ZrZn2, especially given that the 
same electrons are believed to participate both in ferromag- 
netism and SC, this latter scenario had to be taken seriously to 
explain the microscopic coexistence between the two phases. 
In particular, the very large hyperfine magnetic molecular field 
in these materials, measured'-' e.g. with Mossbauer spec- 
troscopy, far exceeds the Pauli limit. This excludes any possi- 
bility of singlet-pairing superconductivity. 

We should note that although the latter statement is true in 
UGe2 and other ferromagnetic superconductors, one may still 
ask whether in principle a singlet-type superconductivity can 



coexist with ferromagnetism. Although some early theoreti- 
cal studies''* indicated that the answer to this question may be 
affirmative provided FM is weak, a more careful analysis con- 
cluded'^ that the coexistence state of spin-singlet pairing and 
ferromagnetism always turns out to be energetically unfavor- 
able against the non-magnetic superconducting state even if a 
finite-momentum pairing (FFLO) state is considered. Later, 
it was proposed'^ that the coexistence of metallic ferromag- 
netism and singlet superconductivity may be realized assum- 
ing that the magnetic instability is due to kinetic exchange. 
However, the coexistence of magnetism and spin-triplet super- 
conductivity appears to be a more promising scenario, since 
the Cooper pairs may use their spin degree of freedom to align 
themselves with the internal magnetic field. 

An experimental fact that is even more striking is that in all 
ferromagnetic superconductors known to date, the SC phase is 
only observed in a small part of the phase diagram otherwise 
occupied by ferromagnetism^^, and it is the region where the 
magnetism appears to be at its weakest that SC sets in - on 
the boundary with paramagnetism when the Curie tempera- 
ture is driven to zero (typically by applying pressure). This 
immediately raises the question of the microscopic origin of 
SC pairing, and whether ferromagnetic spin fluctuations play 
the role of a "glue" for Cooper pair formation very much as 
they do in superfluid ''He. It is equally interesting what role 
the zero-temperature pressure-tuned phase transition plays in 
formation of superconductivity and whether notions involving 
quantum criticality (provided the phase transitions are second 
order) are necessary to explain the phenomenon. 

Although there is no universal answer to this question yet 
and the research efforts, both experimental and theoretical, are 
focused on this issue, it is interesting to note that in both UGe2 
and ZrZn2 the feiTomagnetic phase transition as a function of 
pressure becomes 1st order as the "critical pressure" is ap- 
proached at r = 0. One cannot therefore straightforwardly 
apply a theory of quantum criticality (be it the Hertz-Millis'S 
theory or one of its variations) given the absence of the quan- 
tum critical point as such. It is undeniable however that the 
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point where Curie temperature goes to zero is of crucial im- 
portance to the formation of the SC state. 

Drawing further parallels between triplet-pairing FM super- 
conductors and the superfluid '^He, one may wonder whether 
different symmetries of the SC gap can occur, as is the case in 
the different phasesi^ of "^He. For example, can the gap sym- 
metry with point or line nodes be realized in the ferromagnetic 
superconductors? Very recently, experimental evidence has 
appeared which suggests that the answer is 'yes'. Harada et 
al}^ reported on ^■^Ge nuclear-quadrupole-resonance exper- 
iments performed under pressure, in which the nuclear spin- 
lattice relaxation rate revealed an unconventional nature of su- 
perconductivity implying that the majority spin band in UGe2 
was gapped with line nodes, while the minority spin band re- 
mained gapless at the Fermi level. 

Motivated by this, we present a mean-field model for coex- 
isting ferromagnetism and spin-triplet superconductivity with 
a SC order parameter that displays line nodes in majority- 
spin channel and is gapless for minority spin. We first study 
the interplay between the magnetic and superconducting order 
parameters, and then proceed to make several predictions for 
experimentally relevant quantities: the specific heat capacity. 
Knight shift, and tunneling conductance. Let us briefly sum- 
marize our main results. We find that the low-temperature 
specific heat capacity Cy shows power-lawer behaviour (to 
be contrasted with the conventional exponential decay in the 
s-wave case), and that the gapless minority spins dominate 
the contribution to Cy at low temperatures, giving rise to a 
linear T-dependence. Also, the relative jump in Cy shows 
a strong dependence on the exchange splitting in the system. 
With regard to the Knight shift, we find that it is suppressed 
at T = with increasing exchange splitting of the major- 
ity and minority spin bands when the external field is applied 
perpendicular to the spin of the Cooper pairs in the system. 
In general, however, it depends strongly on the orientation 
of the field with respect to the crystallographic axes of the 
compound, indicative of the triplet pairing in the system. Fi- 
nally, the normalized tunneling conductance spectra show a 
strong directional dependence with respect to the orientation 
of the superconducting order parameter in reciprocal space, 
but change very little upon modifying the exchange splitting 
in the system. Our findings should be useful for comparison 
with experimental studies, and could lead to further insights 
as regards the nature of the superconducting order parameter. 

This paper is organized as follows. We first describe the 
phenomenological framework to be used in this work in Sec. 
im We then present our theoretical model in Sec. |llll and pro- 
vide the results of the self-consistent mean-field treatment at 
both zero and finite temperatures in Sec. |IV A| We then pro- 
ceed to make predictions for experimentally accessible quan- 
tities in Sec. IIVBI using the self-consistently obtained results 
from Sec. |IV A| We discuss our findings in Sec. |Vll and sum- 
marize in Sec. IVIII We will use boldface notation for vectors, 
.?. for operators, .7. for 2x2 matrices, and .?. for 4x4 matrices. 



II. PHENOMENOLOGICAL FRAMEWORK 

The issue of coexisting ferromagnetism and superconduc- 
tivity dates back to half a century ago when the celebrated 
FFLO state was predicted'" as a finite-momentum pairing 
state with real-space structure of the singlet SC order parame- 
ter that may develop under certain conditions close to the crit- 
ical magnetic field Hc2- The conditions for the FFLO state are 
however very different from those observed in ferromagnetic 
superconductors such as UGe2. In particular, as has already 
been emphasised above, the magnetic molecular field felt by 
Cooper pairs inside the ferromagnet is many times larger— 
than the Pauli limiting field necessary to destroy the singlet 
Cooper pairs. We shall therefore concentrate on triplet- type 
superconducting pairing. 

Several remarks are in order We note from the outset that 
the ferromagnetism observed in the Uranium compounds and 
in ZrZn2 is itinerant, Stoner-like in its nature. We shall there- 
fore not discuss the topic of localized magnetic moments that 
would have, among other things, provided a pair-breaking 
mechanism in accord with Abrikosov-Gor'kov theory^' of 
magnetic scattering. Here, we will assume that the same elec- 
trons involved in the spontaneous SU(2) symmetry breaking 
associated with ferromagnetism, also participate in the U(l) 
gauge symmetry breaking that characterizes a superconduc- 
tor 

The idea of triplet pairing occurring between the same elec- 
trons that form the Stoner instability at the border of ferro- 
magnetism goes back to Fay and Appel" (1980) who consid- 
ered exchange of magnetic spin fluctuations as a microscopic 
mechanism for Cooper pairing. More recently, the problem 
has been revisited22i22i2i25i2'' in the light of experimental find- 
ings in UGe2 and other ferromagnetic superconductors. 

In this paper, we shall take a phenomenological approach to 
superconductivity, leaving the intriguing and debated question 
of the microscopic mechanism for Cooper pairing aside. In 
particular, we shall consider systems where superconductivity 
appears at a lower temperature than the temperature at which 
onset of ferromagnetism is found. This is certainly the case 
experimentally and may be simply due to the fact that the en- 
ergy scales for the two phenomena are quite different, with the 
exchange energy naturally being the largest. It may, however, 
also be due to the fact that superconductivity is dependent on 
ferromagnetism for its very existence. Such a suggestion has 
recently been put forth^-. 

A crucial issue to address in this context is whether super- 
conductivity and ferromagnetism are phase-separated (such as 
e.g. solid and liquid phases coexisting at the melting point) 
or not. Fairly strong experimental evidence for non-phase- 
separated coexistence of ferromagnetism and superconductiv- 
ity has recently been presented in UGe2 — . However, even 
if such non-phase-separated coexistence is established, there 
still remains the issue of whether the superconducting order- 
parameters exhibits spatial variations, precisely due to its non- 
phase-separated coexistence with ferromagnetic order One 
obvious candidate for such spatial variations-'"^ is a sponta- 
neously formed Abrikosov vortex lattice, induced by the in- 
ternal magnetization M. As argued in Ref. '35', an impor- 
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tant factor with respect to whether a vortex lattice appears 
or not could be the magnitude of the internal magnetization 
M. Specifically, Ref. |3_6 suggested that vortices may arise 
if 47rM > Hci, where Hd is the lower critical field. It is 
conceivable that a weak ferromagnetic state coexisting with 
superconductivity may give rise to a domain structure, in the 
absence of an external field, that is vortex-free. Therefore, we 
shall consider non-phase-separated coexistence of the FM and 
SC order parameters from here on, as have other studies^^. We 
will also leave the complications arising from the spatial vari- 
ation of the superconducting order parameter originating with 
a putative spontaneously formed Abrikosov vortex lattice in 
the superconducting order parameter for future investigations. 

Spin-triplet superconductors are characterized by a multi- 
component order parameter, which for the simplest case of 
the p-wave may be expressed in terms of three independent 
components of a d-vector: 



^kTT 



-i(A: 



^kTT 



A' 



kTi 



(1) 



Note that dk transforms like a vector under spin rotations. In 
terms of the components of dk, the order parameter itself is a 
2x2 matrix that reads 



Aapik) = (ck,aCk,/3) = [j(dk • Cr)(Jy]a/3, 



(2) 



where cr is the vector of Pauli matrices, and cj^. ^, Ck,a are the 
usual electron creation-annihilation operators for momentum 
k and spin a. 

The superconducting order parameter is characterized as 
unitary if the modulus of the gap is proportional to the unity 
matrix: (A • A^^) cx 1. Written in terms of the vector dk, 
this condition is equivalent to the requirement that (Sk) ~ 0, 
where we have introduced the net magnetic moment (or spin) 
of the Cooper pair 



(Sk) = i(dk X dl) 



(3) 



The unitary triplet state thus has Cooper pairs with zero mag- 
netic moment, whereas the non-unitary state is characterized 
by non-zero value of (Sk) ^ 0. The latter effectively means 
that time-reversal symmetry is spontaneously broken in the 
spin part of the Cooper pairs-^. It is thus intuitively clear that 
having the spin of the Cooper pair aligned with the internal 
magnetic field of the ferromagnet can lower the energy of the 
resulting coexistence state. The above argument that the order 
parameter in the ferromagnetic superconductors must be non- 
unitary has been put forward by Machida and Ohmi^^, and 
others-'*'-**. Distinguishing between unitary and non-unitary 
states in ferromagnetic superconductors is clearly one of the 
primary objectives in terms of identifying the correct SC order 
parameter. To this end, recent studies have focused on calcu- 
lating transport properties of ferromagnetic superconductors 
40.41,42,43.44. 45^ There have also been investigations of identify- 
ing spin-triplet pairing in quasi- ID materials^^i^idM^. 

Finally, we note that inter-subband pairing is expected to 
be strongly suppressed in the presence of the Zeeman splitting 
between the f, I conduction sub-bands. In other words, only 
electrons within the same subband will form Cooper pairs (the 



so-called equal-spin pairing) and we shall set A|| = in 
what follows. Moreover, the requirement of non-unitarity of 
the order parameter then reduces to the requirement that the 
vector dk in Eq. ^ should have two non-zero components, i. 
e. A|| 7^ ^11, which one would expect anyway in the pres- 
ence of the Zeeman splitting between the two spin subbands. 
The spin oftheCooperpair is then (5*2) = i(| A|| p — | A|| p) 
and is aligned along the magnetic field (z being the spin quan- 
tization axis). 



III. THEORY 

We consider a model of a ferromagnetic superconductor de- 
scribed by uniformly coexisting itinerant ferromagnetism and 
non-unitary, spin-triplet superconductivity. We write down 
a weak-coupling mean-field theory Hamiltonian with equal- 
spin pairing Cooper pairs and a finite magnetization along the 
easy-axis similar to the model studied in Refs.i25a26,> namely 

where bi^^ra — (c-ko-Cko-) is the non-zero expectation value of 
the pair of Bloch states. Applying a standard diagonalization 
procedure, we arrive at 

^ = ^o + II^ka7L7ka- 

kcr 

^0=2 X!^^'^'^ ^ ^^'^ ^ A^^^fekffff) H ^ — 

k<T 



, (5) 



where {%a,ltcr} ^^'^ ^^"^ fermion operators and the eigen- 
values read 



^k. = A/C + |Ak..P. 



(6) 



It is implicit in our notation that Ck = Ck — Ep is measured 
from the Fermi level, where £k is the kinetic energy. The free 
energy is obtained through 



(7) 



k(T 



such that the gap equations for the magnetic and supercon- 
ducting order parameters become^^ 



M 



1 y ^tanh(/3£;k./2), 

kcr 



2Eu', 



'-tanh{pEk'a/2). (8) 



Specifically, we now consider a model which should be of rel- 
evance to the ferromagnetic superconductor UGe2, and possi- 
bly also for UCoGe and URhGe. In Ref. |20i, it was argued 
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that the majority spin (spin-up in our notations) fermions were 
gapped and that the order parameter displayed line nodes, 
while the minority (spin-down) fermions remained gapless 
at the Fermi level in the heavy-fermion compound UGe2. 
An obvious mechanism for suppressing the superconducting 
instability in the minority-spin channel as compared to the 
majority-spin channel is the difference in density of states 
(DOS) at the Fermi level. Indeed, from Fig. 1 in Ref. 25| 
(see also Fig. 4 in Ref. |26|) . it is seen that the critical tem- 
perature for pairing in the minority-spin subband, T^, is pre- 
dicted to be much smaller than the critical Tj for the majority- 
spin subband, even for quite weak magnetic exchange split- 
tings. Given the already quite low critical temperature Tc 
that is observed experimentally in ferromagnetic supercon- 
ductors (Tc < 1 K), which we associate with Tj, we there- 
fore conclude that it might indeed be very hard to observe 
experimentally the even smaller gap in the minority-spin sub- 
band. Therefore, it is permissible to only consider pairing in 
the majority-spin channel and neglect a small (if any) pairing 
between minority-spin electrons. In our notation this means 
setting M ^ 0, Akt ^ 0, Akx = 0. 

We stress that the above statement, although intuitively at- 
tractive, may need further justification since we have so far ne- 
glected completely the spin-orbit interaction that is expected 
to be strong in Urainium based compounds, such as UGe2, 
URhGe and UCoGe. The effect of the latter would be to pro- 
vide some effective coupling between majority and minority 
spin subbands and would probably lead to induced SC order 
parameter in minority spin channel. This issue is left for fu- 
ture study^S. 

To model the presence of line nodes in the order parameter, 
we choose 



^kT = Ak^i^ = Aocos6', 



(9) 



where is the normalized Fermi wave-vector, such that 
the gap only depends on the direction of the latter This is 
the weak-coupling approximation. The above gap satisfies 
the correct symmetry requirement dictated by the Pauli prin- 
ciple, namely a sign change under inversion of momentum, 
6 ^ TT — 9. Here, 9 is the azimuthal angle in the xy-plane. 
Our choice of this particular symmetry for the p-wave su- 
perconducting gap is motivated by the experimental results 
of Harada et al^. The cos 0-dependence is also in accord 
with the results of Ref. ^33, which showed that the major- 
ity band at the Fermi level for UGe2 is strongly anisotropic 
with a small dispersion along the A;j,-direction. We consider 
here a situation where the electrons are restricted from mov- 
ing along the z-axis. The motivation for this is that, strictly 
speaking, it seems plausible that uniform coexistence of ferro- 
magnetic and superconducting order should only be realized 
in thin-film structures where the Meissner (diamagnetic) re- 
sponse of the superconductor is suppressed for in-plane mag- 
netic fields. The thin-film structure would then also suppress 
the orbital effect of the field. In a bulk structure, as considered 
in Ref. 14, we expect that a spontaneous vortex lattice should 
be the favored thermodynamical state^**, unless prohibited by 
a possible domain structure. Having said that, we point out 
that there is no firm experimental evidence for the presence 



of such a vortex phase in ferromagnetic superconductors such 
as UGe2 and ZrZn2, and we therefore do not exclude some 
mechanism that would instead stabilise a truly uniform co- 
existence of the SC and FM in these materials. It should be 
mentioned that uniform coexistence of ferromagnetism and 
superconducting order have also been speculated to occur in 
quasi-lD and quasi-2D materials such as RuSr2GdCu208^ 
In our model, the pairing potential may be written as 



V{e,9') = -goose cose', 



(10) 



where g is the weak-coupling constant. Conversion to integral 
equations is accomplished by means of the identity 

^E/(^k.)= / deN^is), (11) 



where N'^{e) is the spin-resolved density of states. In three 
spatial dimensions, this may be calculated from the dispersion 
relation by using the formula 



V 



(27r)3 Jek.=const IVkffko 



(12) 



With the dispersion relation ^kcr = £k ^ alM — Ep, one 
obtains 



^ ^ ^V^'^rn{e + aIM + EF) _ ^^^^ 
In their integral form, Eqs. ([8]l for the order parameters read 



An 



X tanh[/3£;^(e,6i)/2], 

d^de i:"-*.^"^'^ tanh[/3£;| {e,9)/2]. 

(14) 



J -ujo 



For ease of notation, we also define 

if a 



v/e^ + A2 cos2 9. if(T=t 
e if tr =i 



(15) 



For the following treatment, we define M ~ IMj Ep, i.e. the 
exchange energy scaled on the Fermi energy. Moreover, we 
set c = gN {0)/2 to a typical value of 0.2 and a)o = ujq/Ep ~ 
0.01 as the typical spectral width of the bosons responsible 
for the attractive pairing potential. Finally, we define the pa- 
rameter / = IN{Q) as a measure of the magnetic exchange 
coupling. As discussed below, only for I > 1 will a sponta- 
neous magnetization appear in our model, in agreement with 
the Stoner criterion for itinerant ferromagnetism. 
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IV. RESULTS : MEAN-FIELD MODEL FOR 
COEXISTENCE 

A. Zero temperature case 

For zero-temperature, the superconducting gap equation 
reads 

i^JlT r aoae^M^. (16) 

Under the assumption that uq ^ Aq, we obtain that 

— ^ = lnf%^) - - / deicos^OTn I 003 6*1. (17) 
cVl + M ^ ^0 / TT Jo 

which may be solved to yield the zero-temperature gap 



Ao = 2.426cJo exp[-2/(c\/ 1 + M)] . 



(18) 



By inserting Eq. ( fTSl l into the gap equation for the magne- 
tization in Eq. (fT4l l. we have managed to decouple the self- 
consistency equations for M and Aq. Numerical evaluation 
reveals that the gap equation for M is completely unaffected 
by the presence of Aq, which physically means that the mag- 
netization remains unaltered with the onset of superconduc- 
tivity. This is reasonable in a model where the energy scale 
for the onset of magnetism is vastly different from the energy 
scale for superconductivity, such that by the time supercon- 
ductivity sets in, the ordering of the spins essentially exhausts 
the maximum possible magnetisation. 

The dependence of Aq on / is shown in Fig. [T] The gap 
remains constant for / g [0, 1], which is a unitary phase. In 
the unitary phase, there is no reason for the minority spin band 
to remain ungapped when M = 0, and hence we would ex- 
pect two gaps A| = A I of equal magnitude for / < 1. Our 
model of gapping exclusively for the majority spin band is 
therefore justified only for / > 1, which is the regime we 
shall be concerned with throughout this article. The onset 
of a spontaneous magnetization for / > 1 is the well-known 
Stoner criterion for an isotropic electron gas, where the spin 
susceptibility may be written as^ 



1 -ixo{(i,uj)' 



2 



2vF\q\J' 



|q| < 2kF, uj < Ep. 



(19) 



For a parabolic band, the static susceptibility is maximal for 
q = where 



X(q = 0,c^ = 0) = 



iVn 



l-IN{0) 



No 



(20) 



The introduction of a ferromagnetic order is demarcated by 
the divergence of the susceptibility for 1=1, which is pre- 
cisely Stoner's criterion for itinerant ferromagnetism. In the 



absence of superconductivity, the self-consistency equation 
for the magnetization at T = reduces to 



^[(^F + 'jh + r/)3/2 - 2{Ep + ah)^/^], 



(21) 



where t] is an upper energy cut-off determined by the band- 
width and h = IM is the exchange splitting of the majority 
and minority bands. Since the energy scales for the magnetic 
and superconducting order parameter differ so greatly in mag- 
nitude, Eq. (l2n i is an excellent approximation even in the co- 
existent state (we have verified this numerically). 
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FIG. 1 : (Color online) The gap-dependence on the ferromagnetic ex- 
change interaction parameter / = IN{0). The gap remains constant 
for I £ [0, 1], corresponding to a unitary phase. The gap Ao then 
starts growing with increasing / for / > 1.0, announcing the onset 
of a spontaneous magnetization. The analytical formula is based on 
Eq. (till. 



B. Finite temperature case 

The critical temperature for the superconducting order param- 
eter is obtained in the standard way [setting Aq = in Eq. 
(fT6b l to yield 



Tc = 1. 134^^0 exp[-2/(cv/l + M)]. 



(22) 



In Fig. |2l we plot the temperature-dependence of the self- 
consistently obtained solution of Aq and compare it to the 
analytical mean-field temperature dependence 



Ao(r) = Ao(0)tanh[7VT,/r-l]. 



(23) 



The BCS result is 7 = 1.74, but we find a better fit for 
our numerical results using 7 = 1.70. Throughout the rest 
of this paper, we shall therefore make use of Eq. ( |23] ) with 
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^ = 1.70 to model the temperature-dependence of the gap for 
/ = {1-01, 1.02, 1.03}, since the agreement is excellent with 
the full numerical solution. As in the zero-temperature case, 
we find that the gap equations in Eq. (fl4] i may be completely 
decoupled also at finite temperature. We have verified that the 
gap equation for the superconducting order parameter has a 
unique non-trivial solution, which guarantees that the system 
will prefer to be in the coexistent state of ferromagnetism and 
superconductivity. The phase-diagram of the model we are 




FIG. 2: (Color online) Self-consistently obtained solution for the su- 
perconducting gap Ao (red symbols) compared to the analytical ex- 
pression Eq. J23t with 7 — 1.70 (blue lines), modelling a BCS-like 
temperature dependence. 



considering may be obtained numerically and is shown in Fig. 
|3] As seen, a quantum phase transition may occur at / = 1.0, 
separating the 'unitary' superconducting state (see discussion 
in an earlier paragraph) from the ferromagnetic, non-unitary 
superconducting state. The critical temperature for the mag- 
netic order parameter is orders of magnitudes larger than Tc 
for the superconductivity except for very close to / = 1.0. 
The increase in Tc in the non-unitary phase as compared to 
the unitary phase is a result of the increase in density of states 
with magnetization for the majority spin. 

Experimentally, one often maps out the T-p phase diagram, 
where T is temperature and p is pressure. Note that the value 
of / may be controlled experimentally by adjusting the pres- 
sure on the sample. A change in pressure is accompanied by 
a change in the width of the electron bands, and therefore di- 
rectly affects the density of states at the Fermi level: increas- 
ing the pressure on the samples reduces the density of states, 
and hence also the effective coupling constant I.^^. A notable 
feature in the phase diagram for UGe2 as determined experi- 
mentally, is that superconductivity only appears in the ferro- 
magnetic phase, and not in the paramagnetic phase. 
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FIG. 3: (Color online) The phase-diagram of our model in the T-I 
plane. For / > 1.0, a spontaneous magnetization arises and allows 
for the possible uniform coexistence of ferromagnetism and triplet 
superconductivity. Note that decreasing / (going from left to right 
along the x-axis) corresponds to an increasing external pressure p. 
The abbrevations stand for non-unitary (NU), unitary (U), ferromag- 
netic (FM), and paramagnetic (PM). 



V. RESULTS: EXPERIMENTAL PREDICTIONS 

We next proceed to using the self-consistently obtained so- 
lutions from the previous section to make predictions for three 
experimental quantities that are routinely used to study super- 
conducting condensates: specific heat. Knight shift, and tun- 
neling conductance spectra. We first consider the normalized 
heat capacity, which is defined as 



d9de- 



cosh^ 6*) /2] 



dT ) 



(24) 



Since the critical temperature of M is much larger than the 
critical temperature for Ag in our model, we may safely ne- 
glect dM/dT in the low-temperature regime. Consider Fig. 
IHfor a plot of the specific heat capacity using three represen- 
tative values for /. The general trend with increasing / is an 
increase of the jump of Cv at T = Tc- The physical reason for 
this is that the majority spin carriers will dominate the jump 
in specific heat stronger when the exchange splitting between 
the bands increases, which is in agreement with the results of 
Ref. [2^ Analytically, the relative jump in specific heat may 
be expressed as 



Cv 



T=Ta 




(25) 



It depends on the exchange splitting in the superconductor 
since the contribution from the majority spin carriers will tend 
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to dominate the specific heat when h increases. The low- 
temperature scahng with T bears witness of the line nodes 
in the gap, and is to be contrasted with the more rapidly de- 
caying s-wave case. Also note that the minority spin fermions 
are in the normal state and give a significant contribution to 
the specific heat in form of a linear T-dependence at low tem- 
peratures. If both spin species were gapped with line nodes, 
one would expect a -dependence of the low temperature 
specific heat. 

In the experimental study of the heat-capacity in UGe2 con- 
ducted in Ref.[6^ a peak of the heat-capacity associated with 
the superconducting transition was observed in a narrow pres- 
sure region Ap ~ 0.1 GPa around p^. Here, is the pressure 
at which the superconducting transition temperature Tc shows 
a maximum value. Farther away from px, the heat capacity 
anomaly was smeared out. In particular, a substantial resid- 
ual value of Cv /T was observed at T ^ 0. The authors of 
Ref. 60 argued that neither the minority band density of states 
at the Fermi level nor the contribution from a self-induced vor- 
tex state would be appropriate to describe this residual value. 
Instead, it might stem from impurities that induce a finite den- 
sity of states at the Fermi level. For an anisotropic supercon- 
ductor like UGe2, the residual value would be highly sensitive 
to such impurities. It is also clear that the observation of sharp 
peaks, similar to the ones we obtain in Fig. H] depend strongly 
on the applied pressure on the superconductor, and in particu- 
lar to how close it is to p^- 




FIG. 4: (Color online) The left panel shows a plot of the specific 
heat capacity, using self-consistently obtained order parameters, for 
three different values of I. The right panel shows relative jump (su- 
perconducting vs. normal state) of the specific heat at the transition 
temperature as a function of the normalized exchange splitting be- 
tween the spin bands, M. Numerically calculated values are shown 
in red, analytical result [Eq. J25b 1 using 7 « 1.70 are shown in blue. 

We next consider the spin susceptibility, making use of the 
standard formula^' 

X{^.:<^) X! Tr{^(k,itj„)^(k + q,itj + iw„)}, 

(26) 

where Q is the matrix Green's function in particle-hole and 
spin-space, where w„ = 2(n+l)7r//3 are fermionicMatsubara 



frequencies. In the static (uj = 0) and uniform (q = 0) limit, 
Eq. ( l26b reduces to the Knight shift k = x(0, 0). We define 
the normalized Knight shift as 



E 



d6'deiV^(e) 



Ef- 



(27) 



The Knight shift is a measure of the polarizibility of the con- 
duction electrons in the compound, and serves as a highly use- 
ful probe to distinguish between singlet and triplet supercon- 
ductivity. For a singlet superconductor, the total spin S of the 
Cooper pair is zero, and the Knight shift therefore vanishes at 
T = since there are no quasiparticle excitations in the super- 
conductor that may be polarized. The Knight shift vanishes 
regardless of the direction in which the external magnetic is 
applied for a singlet superconductor. For a triplet supercon- 
ductor, this is quite different. The Knight shift now may be 
anisotropic in terms of the direction in which the magnetic 
field is applied. By means of the dk-vector formalism [see 
Eq. ([1])], one may infer that the Knight shift is unaltered even 
for T < Tc when dk -L H, but is altered according to Eq. [27] 
when dk| |H. This is valid as long as the dk remains 'pinned' 
in the material due to e.g. spin-orbit coupling, and hence does 
not rotate with H. Otherwise, the Knight shift would remain 
unaltered in any direction. Therefore, an anisotropic Knight 
shift is a strong signature of a vector character of the super- 
conducting order parameter, and hence of a spin-triplet super- 
conducting state. 

In Fig. |5] we plot the Knight shift for several values of /. 
It is interesting to note that k(0) is reduced with increasing /. 
Physically, this may be understood by realizing that the den- 
sity of states of ungapped minority spins at the Fermi level 
decreases as the exchange splitting between the majority- and 
minority bands increases. This results directly in a lower 
amount of polarizable quasiparticles, and hence the Knight 
shift becomes suppressed. For a fully polarized ferromagnet 
(half-metal), the Knight shift would therefore be identical to 
an .s-wave singlet superconductor for an applied field satisfy- 
ing H 11 dk. This fact emphasizes the importance of measur- 
ing the spin susceptibility along several directions to identify 
the proper spin-symmetry of the superconductor 

As a final experimental probe for the interplay between fer- 
romagnetism and superconductivity, we employ a Blonder- 
Tinkham-Klapwijk formalism.^'^ to calculate the tunneling be- 
tween a normal metal and a ferromagnetic superconductor 
in the clean limit, using the self-consistently obtained values 
of the order parameters in the problem. From the results of 
Ref. [43, we find that the normalized tunneling conductance 
may be written as 



Go 



E 



7r/2 



Tr/2 



decose[l + \r^{eV,t 



„N 



(28) 



where Gq is the normal-state conductance. Above, r^{eV, 9) 
and (eV, 0) designate the Andreev- and normal-reflection 
coefficient, respectively, and read 
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FIG. 5: (Color online) Knight shift for a ferromagnetic supercon- 
ductor, using self-consistently obtained order parameters, for three 
different values of I. The field is here applied H || dk- 



u.{0uh.m-)m\^ 



4fcf cos0q- cos0^v„{0-,+ )u,{0^_)rA&U) 



(29) 



We have defined Z — 2mVb/fcF as a measure of the barrier 
strength, where m is the quasiparticle mass, Vq is the scatter- 
ing strength of the barrier, and kp is the Fermi momentum. 
Moreover, is the angle of incidence of incoming electrons 
from the normal side and we have implicitly incorporated con- 
servation of group velocity and conservation of momentum 
parallel to the barrier, i.e. hp sin = q'^ sm0'^. Finally, we 
have introduced 

T"^^ q" cos 0"^ ±kF cos 0±ikFZ (30) 

and 7,(0) = A,{0)/\A,{0)\, 0-^ - 0-, 0^_ = tt - 0- . In 
the quasiclassical approximation Ep ^ (Ao,e), the wave- 
vectors read 

kp = ^2mEF, q" = ^ImiEp + alM) (31) 
while the spin-generalized coherence factors are 

u.(0%^) = -i={i + ^\-{\A„{0i^)\iEYYi\ 

vA9:±) - ^{1 - .^l-(|A,(0J±)|/i?)2}V2. (32) 

In Fig. |6] we plot the conductance spectra of a nor- 
mal/ferromagnetic superconductor junction. By writing the 
gap as A = Ao cos(0 — a), we allow for an arbitrary orienta- 
tion of the gap with respect to the crystallographic axes. The 



features seen in the conductance spectra are qualitatively dif- 
ferent for a — and a = 7r/2. In the first case, the electron- 
and hole-like quasiparticles entering the superconductor ex- 
perience a constructive phase-interference which gives rise to 
the formation of a zero-energy state that is bound to the sur- 
face of the superconductor. The resonance condition for the 
formation of such zero-energy states is A(0) = — A(7r — 0),— 
and the bound states are manifested as a giant peak in the 
zero-bias conductance^"*. Note that such states exist even if 
the spatial depletion of the superconducting order parameter 
is not taken into account, which may be shown analytically^-. 
Taking into account the reduction the gap experiences close to 
the interface compared to its bulk value, is known to yield 
the same qualitative features as the usual step-function ap- 
proximation, with the exception of additional, smaller peaks 
at finite bias voltages due to non-zero bound states^^. From 
Fig. |6] we see that the effect of increasing the exchange field 
amounts to sharper features in the conductance spectra. With 
increasing /, the zero-bias conductance peak becomes larger 
for a = 0, while the dip structure for a — tt/2 becomes 
more pronounced. Physically, this may be understood by the 
increased contribution from majority spin carriers. The con- 
tribution from the minority spin carriers is constant for the en- 
tire low-energy regime, and leads to less pronounced features 
in the conductance. The effect of the barrier strength Z is seen 
in the left column of Fig. |6] For a — 0, increasing Z leads to 
a higher peak at zero bias, while increasing Z suppresses the 
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FIG. 6: (Color online) Plot of the tunneling conductance of a normal/ferromagnetic superconductor junction for a = Q and a — 7r/2, using 
self-consistently obtained solutions at T = 0. In the left column, we fix the tunneling barrier strength Z = 2mVo/kF = 3 and plot the 
conductance for several values of the Stoner interaction /. In the right column, we fix J = 1.01 and plot the conductance for several values of 
Z. 



conductance for a — tt/2. 

It is also worth emphasizing the relation between the tun- 
neling conductance and the bulk DOS of the superconductor 
As is well-known, the conductance of a normal/s-wave super- 
conductor junction in the tunneling limit approaches the DOS 
of the bulk superconductor^^. The same argument is valid for 
a c?a,2_y2-wave superconductor^"*. One might be tempted to 
conclude that the tunneling conductance will always approach 
the bulk DOS of the superconductor in the strong barrier limit 
as long as there is no formation of zero-energy states. How- 
ever, closer examination reveals that this is not necessarily so. 

To illustrate this, we draw upon some results obtained in 
Ref.Ugj. In general, the conductance of an N/S junction in the 
tunneling limit may be written as 



PS = Po, where 



GieV) 



^6* AT (TAT cos 9NPs{eV) 



'^^N(7N COS e 



(33) 



N 



where crjv is the normal-state conductance for a given angle 
of incidence 9^ and ps is the surface DOS for the super- 
conductor. In the absence of zero-energy states, the surface 
DOS coincides with the bulk DOS of the superconductor, i.e. 



PoieV) = 



/ d6'7vRe<^ - 



eV 



-V2 ^ ^eV^-\A{9M)? 



}. (34) 



An important consequence of the above equation is that the 
tunneling conductance may be interpreted as the expectation 
value of ps with a weighting factor ctat cos 6'jv. 

Let us now compare three different superconducting sym- 
metries to illustrate the relation between the conductance and 
the DOS. We consider an s-wave, (i^2_j^2-wave, and p^-wave 
symmetry, none of which feature zero-energy surface states 
(Fig. [8]l. Naively, one might therefore expect that the con- 
ductance should converge towards po in the tunneling limit. 
However, it turns out that the weighting factor gn cos 9n, 
which is peaked around 9n = 0, plays a major role in this 
scenario. In Fig. [T] we plot both the tunneling conduc- 
tance G{eV)/GQ and the bulk DOS po for these three sym- 
metries and fix Z — 20. We regain the well-known results 
that G{eV)/GQ po for large Z in the s-wave and d,j.2_y2- 
case. However, the conductance and DOS differ in the py- 
wave case. 

The reason for the deviation between G/Gq and po in the 
Pj,-wave case may be understood by consulting Fig. [8] As 
seen, the weighting factor is peaked around normal incidence 
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FIG. 7: (Color online) Plot of the normalized conductance G/Go 
and bulk DOS po for three different symmetries of the superconduct- 
ing state in the tunneling limit. Only in the p^-wave case is there a 
difference between these two quantities. 
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FIG. 8: (Color online) Illustration of the different symmetry states 
considered here and a qualitative sketch of the weighting factor 

CTJV COS 9n- 



f?Ar — 0. In the s-wave and d,j.2_y2-case, the gap magnitude 
is maximal at 6'jv = and replacing the weighting factor in 
Eq. ( |33T l with unity has little or no consequence. The situation 
is dramatically different in the pj,-wave case. Now, the gap 
magnitude is actually zero for normal incidence, and it is pre- 
cisely this contribution that will dominate the integration over 
angles in Eq. ( [33T l. Therefore, replacing the weighting factor 
with unity, in order to obtain the DOS, has a non-trivial con- 
sequence in the pj,-wave case. This analysis illustrates how 
the conductance and bulk DOS in the absence of zero-energy 
states are not always the same in the tunneling limit. Note that 
the orientation of the interface with respect to the symmetry 
of the order parameter is crucial with regard to the measured 
conductance spectra and the surface DOS. For instance, even 
at Of = 7r/4 there is an appearance of a large zero-bias conduc- 
tance peak for the p-wave pairing considered here, although 
the gap orientation does not satisfy the condition for perfect 



formation of zero-energy states. 



VI. DISCUSSION 

We have discussed a mean-field model where itinerant fer- 
romagnetism coexists with non-unitary, triplet superconduc- 
tivity, with a gap that contains line nodes. The precise sym- 
metry of the order parameter in the ferromagnetic super- 
conductors UGe2, URhGe, UCoGe is still under debate, al- 
though most experimental findings and theoretical consider- 
ations strongly point towards the realization of a triplet su- 
perconducting order parameter It is plausible that such a su- 
perconducting order parameter is non-unitary, thus breaking 
time-reversal symmetry in the spin channel of the Cooper pair 

The orbital symmetry of the superconducting order parame- 
ter in ferromagnetic superconductors is a more subtle issue. In 
Ref. a mean-field model for isotropic, chiral p-wave gaps 
in a background of itinerant ferromagnetism was constructed. 
In that work, pairing was assumed to occur both for majority- 
and minority-spins, resulting in for instance a double-jump 
structure in the specific heat capacity. An isotropic, chiral 
p-wave order parameter has a constant magnitude, which is 
favorable in terms of maximizing the condensation energy 
gained in the superconducting state. Assuming an isotropic 
density of states at the Fermi level and a separable pairing po- 
tential of the form Vkk' = ^ gAkAks the condensation energy 
gained at T = in the superconducting state reads 



E 



iV(0)A2 



(|Ak|^ 



(35) 



where Aq is the maximum value of the gap and (. . .) denotes 
the angular average over the Fermi surface. This clearly shows 
the advantage of an isotropic gap |Ak| = 1. The general prin- 
ciple is well-known: the system prefers to have the Fermi sur- 
face as gapped as possible. However, factors such as spin- 
orbit pinning energy and lattice structure may conspire to pre- 
vent a fully isotropic gap. We also note that in our model, the 
ferromagnetic ordering enters at a much higher temperature 
than the superconducting order unless / is very close to unity. 
This is consistent''^ with the experimental findings for the ra- 
tio between the critical temperatures for feiTomagnetic and su- 
perconducting order, T™ /Tf~^, except for UCoGe where the 
ratio is ~ ?>A 

The experiments performed so far are indicative of a single 
gap, or at least a strongly suppressed second gap, in the fer- 
romagnetic superconductors. For instance, no double-jump 
features have been observed in the specific heat capacity' for 
UGe2. This warrants the investigation of a single-gap model, 
possibly with line nodes as suggested by Harada et a/. 3 The- 
oretically, the absence of the SC gap in the minority spin 
subband can be justified by considering the effect of Zeeman 
splitting on the electronic density of states (see discussion in 
SecHniand Ref. 25). In general, it should be possible to dis- 
cern the presence of two gaps by analyzing specific heat or 
point-contact spectroscopy measurements, unless one of the 
gaps is very small. 
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Apart from this, another possible scenario, specific to 
UGe2, can be invoked to explain the observed gapless be- 
haviour in the minority spin subband. This is the meta- 
magnetic transition that occurs inside the FM phase of UGe2 
and separates the two ferromagnetic phases with different val- 
^g j60,6i magnetization M. The reason this meta-magnetic 
transition in UGe2 is of great importance is because the spe- 
cific heat measurements clearly indicate*" that the maximum 
of superconducting Tc occurs not at the FM to PM transition, 
but at some lower pressure « 12 kbar that coincides pre- 
cisely with the meta-magnetic transition-"'*'. 

One can think of this transition as a point where the value 
of low-temperature magnetization Af sustains a jump. While 
the microscopic origin of this transition is not known, an idea 
has been put forward*^ that it may be due to a sharp change 
in the density of states (DOS) due to the existence of a double 
peak in its structure close to the Fermi level. What happens 
according to this scenario is that applied pressure makes the 
Fermi level "sweep through" the double-peak structure in the 
DOS, thereby sharply increasing the density of states in the 
majority spin channel. It follows from a simple Stoner insta- 
bility argument that such an increase in the DOS would lead 
to a larger value of effective interaction / = IN{Q) and thus 
higher magnetization M. But this also means that the ratio of 
the DOS in the two spin channels, iV|/A^|, sharply increases 
at the meta-magnetic transition. It follows from Eqs. ( [TSlfTSl 
|22] | that the ratio between the SC gaps in the two spin sub- 
bands 

At tJ cxpi-l/gN^) ^■'"^ 

thus becomes very small, justifying the assumption A| = 
made in this work. 

We note in passing that from an experimental point of view, 
a complication with UGe2 is that the superconductivity does 
not appear at ambient pressure, in contrast to URhGe and 
UCoGe. The necessity of considerable pressure restricts the 
use of certain experimental techniques, and this is clearly 
a challenge in terms of measuring for instance conductance 
spectra of UGe2 . Another experimental quantity which would 
be of high interest to obtain from for instance ab initio cal- 
culations, is the thermal expansion coefficient, which may be 
directly probed in high-pressure experiments^ ^ 

We also underline that in our model the magnetism is as- 
sumed to coexist uniformly with superconductivity. Depend- 
ing on the geometry of the sample, it is likely that the intrinsic 



magnetization gives rise to a self-induced vortex phase. In a 
thin-film structure where the thickness t is smaller than the 
vortex radius A, we expect that ferromagnetism and supercon- 
ductivity may be realized in a vortex-free phase, similarly to 
a thin-film s-wave superconductor in the presence of an in- 
plane magnetic field. Further refinements leading to a more 
realistic model of a ferromagnetic superconductor should in- 
clude the presence of spin-orbit coupling, which inevitably is 
present in heavy-fermion superconductors, in addition to the 
presence of vortices. Nevertheless, we believe that our model 
should capture important qualitative features of how the in- 
terplay between ferromagnetism and superconductivity may 
be manifested in experimentally accessible quantities. In par- 
ticular, experiments on transport properties of ferromagnetic 
superconductors, such as the Josephson current and point- 
contact spectroscopy, would be of high interest to farther illu- 
cidate the pairing symmetry reaUzed in ferromagnetic super- 
conductors. 



VII. SUMMARY 

In conclusion, we have constructed a mean-field theory of 
triplet superconductivity in the background of itinerant ferro- 
magnetism, where the superconducting order parameter con- 
tains line nodes and the minority spin band remains ungapped 
at the Fermi level. We have solved the self-consistency equa- 
tions for the order parameters in the problem, and find that 
ferromagnetism enhances superconductivity, while the ferro- 
magnetism itself is virtually unaffected by the presence of su- 
perconductivity. We have made several predictions for exper- 
imentally accessible quantities: heat capacity. Knight shift, 
and tunneling conductance spectra. Our results may be helpful 
in the interpretation of experimental data, and could provide 
tools concerning the issue of identifying the pairing symmetry 
of ferromagnetic superconductors. 
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